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Abstract 

We analyze the pattern of the onset of complete absorption (the black limit) in the high 
energy hadron-hadron collisions. The black limit arises due to the hard and soft interaction 
dynamics as a function of the impact parameters b. Both hard and soft mechanisms provide 
universal dependence of the partial amplitude of the high energy elastic hadron- hadron scatter- 
ing on the impact parameter b and result in the radius of interaction proportional to ln(s/so). 
We find that with increase of the collision energies hard interactions lead to a faster increase 
of the impact parameter range where the partial wave amplitudes are approaching the unitarity 
limit. Consequently, we argue that at super high energies when the radius of hadronic interac- 
tions significantly exceeds static radii of the interacting hadrons(nuclei) the ratio of total cross 
sections of nucleon-nucleon, meson-nucleon, hadron-nucleus, nucleus-nucleus collisions be- 
comes equal to one. The same universality is also expected for the structure functions of nuclei: 
F2a{x, Q 2 )/F2n(x, Q 2 ) — ► 1, at very small x, and for the ratio ct 7 a/o" 7P at superhigh ener- 
gies. We analyze how accounting for the energy dependence of the interaction radii changes the 
geometry of hadron-nucleus and nucleus-nucleus collisions, the energy dependence of total, 
absorption and inelastic cross sections, the distribution over the number of wounded nucle- 
ons in proton-nucleus collisions and find that these effects are noticeable already for the LHC 
energies and even more so close to the Greisen-Zatsepin-Kuzmin limit. 



1 Introduction 

High energy behavior of the hadron-hadron interaction is a subject of the extensive theoretical 
studies [1J for almost fifty years. Still many long standing questions such as the universality of the 
hadronic and hadron-nucleus total cross sections [2| and the relative importance of hard and strong 
interactions as a function of the impact parameter remain open. 
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V.Gribov [3| demonstrated long ago that the t channel unitarity of the S-matrix combined with 
analytic properties of the scattering amplitude in momentum transfer space is inconsistent with the 
radius of a hadron being independent of the incident energy. An assumption that the single pole 
singularity in the angular momentum plane (a Pomeron pole) determines high-energy behavior of 
the elastic amplitude |@J 15) predicts an increase of the interaction radius with increase of energy. 
This phenomenon leads to increase of the slope parameter B pp in the t-dependence of the pp elastic 
cross section which is by now well established experimentally -see [6| and references therein. 
Estimates within the Pomeron pole model show that the radius of interaction becomes comparable 
with the mean internucleon distances in nuclei already at the LHC energies and even exceeds these 
distances at energies corresponding to the Greisen-Zatsepin-Kuzmin (GZK) limit Q. This leads 
to a gradual change of geometry in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions. 
Also, it results in a gross difference between structure of the final states in central collisions (for 
example triggered by x T = Ip^j^fs > 0.01 dijet production) and peripheral collisions at the 
Tevatron energies and above 1 8 1 . 

Another fundamental property of the strong interaction is a rather rapid increase of the total 
hadron-hadron cross section with energy. Such a behavior is well described, for example, in the 
Donnachie-Landshoff model [9| where the intercept of the Pomeron trajectory is op(0) = 1.08. 
However, the hypothesis of the Pomeron exchange implies the existence of the Pomeron branch 
points in the plane of angular momentum which compete with Pomeron exchange ifTUIl . Measure- 
ments of the total inelastic diffraction cross section in pp scattering ifTTI discovered the important 
role of shadowing effects which can be described as due to the presence of the Pomeron branch 
points (for recent discussions of these issues see Refs. El). Account of the Pomeron branch 
points in the model with the value of ap(0) — 1 > predicted blackening of interaction at the 
range of impact parameters b oc ln(s/s ) (see [T7~3l and references therein). This corresponds to 
the Froissart-type behavior of the total cross section in the s —> oo limit [14|. Experimentally, 
at energies of the Tevatron partial amplitude for pp scattering at impact parameter b = reaches 
a value ~ 0.95. Thus both theory and experimental data indicate that the blackening of the soft 
interaction at particular fixed impact parameters, b, until the unitarity limit is reached is one of the 
distinctive features of soft hadron-hadron dynamics at high energies. 

On the other hand theoretical study of DIS in perturbative QCD (pQCD) lfT?ll and experimental 
investigations at HERA [16] found that the total cross sections of DIS, being small, undergo a 
significantly faster increase with energy than the cross sections of soft hadron-hadron interactions. 
The significant cross section of diffraction observed in DIS at HERA |[T6ll gives convincing evi- 
dence for the blackening of hard QCD interactions at small x ifPTl . In proton-proton collisions at 
collider energies, hard parton interactions are concentrated at significantly smaller impact param- 
eters than generic inelastic interactions 1 8 1 . However, similar to the DIS case, these hard parton 
interactions rapidly increase with energy leading to blackening of interactions in the wide range 
of central impact parameters. Thus it is important to analyze the relative importance of soft and 
hard QCD interactions in the onset of the black body regime of high energy hadron collisions as a 
function of the impact parameter. 

In section 2 we analyze the proton-proton elastic scattering amplitudes as a function of impact 
parameter b. We argue that the black regime originates from hard interactions of leading partons 
in the nucleons with the small x gluon fields. We use this to estimate the rate of increase of the 
region in impact parameter space where the interaction is black and we find that it linearly increases 
with lns/so. We also give an alternative estimate based on the extrapolation of the soft Pomeron 
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dynamics which dominates at large b to the region of smaller b and find that two estimates give 
similar results at collider energies. We use these observations to argue in section 3 that a universal 
pattern of blackening of the interaction at fixed impact parameters leads to the universality of cross 
sections in the limit of superhigh energies: 

crtat(hi,h 2 )/crtot(pp) 1- (1) 

Here hi can be a hadron (nucleus). The same result should be valid for the structure functions of 
nuclei but at extremely small x: 

F 2A (x,Q 2 )/F 2N (x,Q 2 )^l. (2) 

In the string models, similar universality arises as the universality of the term in the cross section 
which increases linearly with energy - the "Pomeron" (loop contribution) related to the universality 
of the graviton interaction [fl~8l . 

It is worth mentioning here that universality of cross sections of hadron-hadron, hadron-nucleus 
collisions at superhigh energies has been discussed long ago. V.Gribov suggested universality of 
cross sections of hadron-hadron, hadron-nucleus interactions [19| within the assumption that such 
cross sections should become constant at infinite energies. Later on, the universality of the Froissart 
limit for total cross sections of hadron-hadron, hadron-nucleus and nucleus-nucleus collisions was 
suggested within a particular generalized eikonal model of the supercritical Pomeron ll20l . This 
model gives prescriptions for summing the unstable (divergent) series of terms oc s n( - aF ~ 1 ^ due to 
multiPomeron exchanges which rapidly increase with energy because of the intercept of the bare 
supercritical Pomeron «jp(0) > 1. There exist also a number of eikonal models which combine 
elements of soft and hard dynamics for all impact parameters, see lETl and reference therein. 
However, universality and the increase of the radius of the hard interactions were not discussed in 
these models. Note also, that the eikonal models neglect contribution of the enhanced Pomeron 
diagrams which is rapidly increases with energy and at superhigh energies becomes comparable 
with eikonal diagrams ll2"2l . 

Our approach assumes dominance of the Donnachie-Landshoff soft Pomeron exchange in pe- 
ripheral collisions only. For the scattering at central impact parameters, where according to pre- 
QCD Reggeon Calculus strong interaction between the Pomerons is expected, our approach ac- 
counts for the large cross section of hard processes due to formation of large gluon densities with 
p± « 2GeV and the consequent disappearance of soft elastic and diffractive processes at energies 
of LHC and above 1 . Such nontrivial interplay of hard and soft dynamics is absent in preQCD mul- 
tiPomeron exchange models. Note here, that importance of hardinteractions in the the Froissart 
limit in QCD follows from requirement of the self consistency ll23l . 

We also analyze how the increase of the radius and strength of the nucleon-nucleon interaction 
influences cross sections of hadron-nucleus and nucleus-nucleus collisions at achievable energies. 
To visualize the role of these effects we use the formulae of the Glauber- Gribov model [24 1 2 . We 

'Disappearance of diffractive processes for the scattering at central impact parameters has been observed recently 
at FNAL. Remember also that in the black limit elastic scattering occur at peripheral impact parameters only. In the 
realistic situation inelastic diffraction arises from the scattering at peripheral impact parameters where the interaction 
is grey. 

2 V.Gribov has demonstrated that though the set of diagrams which contributes at high energies is different from that 
accounted in the original Glauber approach, the answer differs only due to the contribution of the inelastic diffraction 
- the inelastic shadowing corrections. Relative contribution of these corrections is rather modest and decreases at very 
high energies due to suppression of the inelastic diffraction. 
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calculate the total and absorption cross sections for pA collisions. The noticeable effects due to 
the blackening of the nucleon-nucleon interaction can already be observed at the LHC energies 
and above. The onset of black body limit leads to gradual weakening of the A-dependence of 
cross sections which ultimately results in the A-independent cross sections corresponding to the 
universality regime of Eq. [U 

We evaluated the effective energy dependent radius of a nucleus and estimated total cross sec- 
tions for heavy ion collisions using the popular Bradt-Peters expression [25|. We compare this re- 
sult to that obtained with more refined Glauber- like model of nucleus-nucleus interaction ll26ll27l . 
Also we calculate the distribution over the number of inelastic collisions in nucleon-nucleus inter- 
action ll28l and find that accounting for the energy dependent radius of a nucleon-nucleon interac- 
tion leads to a significant change of distribution. Analysis of the possible role of this effect in the 
interpretation of the heavy ion collisions data in terms of the wounded nucleons already at energies 
of RHIC and above is beyond the scope of this paper. 

2 Partial wave amplitudes for hadron-hadron collisions at ul- 
trahigh energies 

In this section we will use the properties of impact parameter representation of the elastic scattering 
amplitude, 

T(b, s) = ^^ J d 2 qexp[-iq- b]/jvjv(q, s), (3) 

to discuss basic features of the superhigh energy NN interactions. Here s is the invariant energy 
of NN scattering, and q is the two-dimensional transverse momentum. 



2.1 Small impact parameter behavior of the hadron collision amplitude 

There are several generic features of T(b, s) at small impact parameters b which can be derived 
from unitarity of S matrix, from the current understanding of the spatial structure of the fast nu- 
cleon and dynamics of hard interactions. Indeed, let us consider nucleon-nucleon scattering at 
small impact parameters at large s. An analysis, performed in Ref.[8|, demonstrates that in this 
case the average transverse momenta of partons from nucleon become large after partons pass 
through the low x gluon fields of another nucleon. For example, at Tevatron energies quarks with 
x > 0.2 get average transverse momenta > lGeV/c. If a leading quark gets a transverse mo- 
mentum pj_, the probability for the nucleon to remain intact is roughly given by the square of the 
nucleon form factor F^^p 2 ^). Since F^(p 2 ± ) < 0.1 for p± > lGeV/c, the probability of survival 
averaged over p± should be at most 1/2, provided average p± > lGeV/ c. Since there are six lead- 
ing quarks (plus a number of leading gluons) the survival probability for two nucleons on small 
impact parameter, | 1 — T(b, s) | 2 , should go as a high power of the survival probability for the case 
of one parton removal, rs (1/2) 6 . Consequently, |1— T(b ~ 0, s)| 2 is close to already at Tevatron 
energies. 

Hence we conclude that 

r(b ~ 0, y/s > 2 TeV) « 1. (4) 

This evaluation is in a good agreement with the Tevatron data on elastic pp collisions and with 
their extrapolations to the LHC energies. 
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Since the probability of hard inelastic interactions at fixed impact parameter increases with 
energy at least as the gluon density at small x, xGn(x, Q 2 ) oc x~ n '\ the increase of this probability 
should be proportional to s Hh with > 0.2. The HERA data suggest that taming of the interaction 
of small dipoles starts only when the probability of inelastic interaction becomes large enough 
(> 1/2))- However for such probabilities of single parton interactions the mulitparton dynamics 
insures the overall interaction to be practically black. Hence, the multiparton interactions will 
ensure a rapid (power law) onset of the regime of black interactions. 

The analysis of the HERA data suggests |29j| that the transverse distribution of the hard 
partons (at the resolution scale p t ) in nucleons can be described as oc exp(— rrih(x)b). Since the 
interaction amplitude of the hard high energy interaction is oc s nh we find that the range of b < hp 
where the interaction is completely absorptive should depend on energy as 

rih In — 

b F « -H-jf, (5) 
m h (x) 

which corresponds to the Froissart limiting behavior. Note that Eq|5] is obtained in the limit of 
sufficiently high energies when rrih(x) for x resolved at the corresponding energy (x ~ Ap 2 / 1 s 
where p t is hard scale) is much smaller than that for the fast partons. 

Actually, a few uncertainties limit our accurate knowledge of the approaching to the black 
body limit due to the mechanism of hard interactions. First at all, there are the uncertainties in 
the x dependence of the gluon densities at very small x where one needs to take into account 
both ln(xo/x), ln(Q 2 /Qo) effects (for the recent discussion see (30)). Also one has to account 
for the increase of the transverse spreading of the gluon distribution with decrease of x. The 
small x evolution is likely to lead to increase of n h at very small x and sufficient virtualities. The 
neglected smearing of the fast partons distribution leads to a decrease of effective mass parameter 
at preasymptotic energies and to a somewhat faster increase of bp with energy. Hence, both effects 
are likely increase the rate of the change of b F at extremely small x. 

The value of m h and rate of its decreasing with energy for x < 1CT 4 can be estimated based 
on the extrapolation of extracted from the HERA data for J/ ip photoproduction which cover 
x > 1CT 4 and correspond to m h (x = 1CT 4 ) ~ 0.75GeV: 

m h (s) = 0.75 [l- 0.027 m(— )]. (6) 

St 

The expected limiting value of is 2m^. Hence, with the reasonable value rih = 0.25 we find 
that the true asymptotic is likely to be reached at fantastically high energies s « 10 22 GeV 2 . At 
these energies m h will be the same for any colliding hadrons build of light quarks. Note that in 
pQCD n h is determined by gluon distribution and, hence, will be also universal, the same for any 
colliding particles. In the case of hadrons with hidden heavy flavor, the onset of universality regime 
requires significantly larger energies. 



2.2 Large impact parameter behavior of the hadron-hadron amplitude 

To determine the behavior of the amplitude at large b we can use arguments based on soft physics. 
At large b only single Pomeron interactions are possible as all multiPomeron interactions have 
much smaller radius (y/2 times smaller for the double Pomeron exchange, etc). Hence we can 
use here information from the analysis of the data on pp scattering at collider energies. Since 
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we are interested in the large b behavior of the amplitude which is determined by the properties 
of the scattering amplitude at small t, we can use the simplest exponential parameterization of 
the t dependence of elastic scattering, leading to the proton-proton amplitude due to the Pomeron 
exchange : 

r ^ b ' s ) = a p ^ ex P [-b 2 /2Bjp(s)} . (7) 
AtcBp(s) 

We take the amplitude to be imaginary at high energies (the ratio of real part of amplitude to the 
imaginary one at high energies we consider here is practically constant and small k « 0.1 -j- 0.13). 
The parameter, 

B P = B 0; p + 2apln(s/s ) (8) 

is the slope of the t dependence of the Pomeron exchange contribution into the amplitude of elastic 
hadron-hadron collision OlEl. The discovery of the diffraction cone shrinkage with increase of 
the energy in the elastic pp collisions confirmed experimentally the energy dependence of Bp 
with the values of parameters B 0t p »i 8.5Gel A ~ 2 ,So = 1 GeV 2 and a'p = 0.25 GeV~ 2 . The 
Pomeron exchange hypothesis predicts Bp/ Bq p « 1.5 at the RHIC energies, Bp/B p « 2 in 
the kinematics of the LHC and Bp/B p rs 2.7 when energy is close to the GZK limit. Thus 
separate analysis of peripheral and central collisions becomes appropriate at energies of LHC and 
above. While the Pomeron model predicts that at extremely high energies the partial amplitude at 
large impact parameters is determined by the Pomeron exchange [31 1, the situation is much more 
complicated at smaller b because of shadowing effects (Pomeron branch points). This will lead 
to renormalization of the Pomeron parameters and to blackening of the interaction if atp(o) > 1 
cf. lfT3ll . So, the calculation of the partial wave amplitude at non-peripheral b is technically rather 
cumbersome. Instead, we assume (as discussed in the Introduction) that partial amplitudes for 
b < bp are Imf(b < bp, s) = 1 and at b > bp are given by the dominant single Pomeron 
exchange 3 . Then, the minimal estimate of bp can be obtained by requiring the continuity in the 
matching of these two regimes: 

pp i \ 

r P (b F ,s) = £f^yex P [-b F /2Bjp(s)] = 1. 

For crfot( s ) we can use the Landshoff-Donnachie parameterization of the Pomeron contribution 
<Jto t (s, s ) = c [^] ap ^ 1 which provides a good description of the data in the region between ISR 
and Tevatron energies with parameter ctjp(0) — 1 = 0.0808. Hence, we get the following energy 
dependence from the condition that at the energies of Tevatron, sp, the partial wave of the NN 
amplitude at impact parameter b = becomes black (T(b = 0, sp) = 0.9 0.95 ~ 1): 

b^ w (a P (0) - 1) Hs/s T )Bp. (9) 

At s — + oo the parameter bp is universal-the same for any colliding particles: 

bp^2a'p{ap{0) - l)ln 2 (s/s T ). 

Another potentially important soft contribution to the partial wave amplitude may arise from 
the hadron scattering off meson tails. Using the dispersion representation of the amplitude over 

3 MultiPomeron cuts are decreasing with increase of b much faster than Pomeron exchange 1 3 1 1. For example, if 
rV(b) cx exp(— ah 2 ), n-Pomeron cut is oc exp(— anh 2 ). 
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momentum transfer, t, it is easy to obtain ll32l : 

r(6) = cs^^expl-fib}. 

The natural expectation given the fast decrease of the amplitude with increase of b is that /x is the 
minimal mass permitted in the channel with vacuum quantum numbers i.e /j = 2m,,-. Thus at 
energies when partial waves with fixed impact parameter b become equal one we obtain: 

b F « l//i(a(/i 2 ) - l)ln(s/s ). 

This value is smaller than that provided by the Pomeron exchange model, 

b F = ^B P (a(0)-l)Hs/sT), 

at achievable energies but exceeds it at asymptotic energies. The leading term in bp is the same 
if the partial amplitude is less than one. It is easy to check that the combination of two types of 
the b dependencies discussed in this subsection does not change the conclusion concerning the 
universal value for the leading term in bp at superhigh energies. 

2.3 Matching of small b and large b behavior 

We demonstrated above that the elastic amplitude is well constrained both at small and large impact 
parameters. To build a complete description we need to determine at what bp the two regimes 
match. Both soft and hard approximations give practically the same result for bp at energies in the 
range from the LHC up to the GZK limit - see FigUJ This is due to the presence of a large constant 
term in Bp. Hence we find that there is very good consistency between the logic of matching 
starting from small b and the logic of matching from large b. It means that there is a smooth 
transition from hard regime at small b < bp to soft regime at higher b > bp for the whole range 
of energies which maybe probed experimentally at colliders and in cosmic ray interactions near 
the GZK cutoff. At the same time we find that the asymptotic rate of the increase of bp is a factor 
of two larger in the hard matching approximation. This indicates that at very high energies, the 
hard mechanism component of the black interactions gives a dominant contribution to the cross 
sections. Also, both hard and soft dynamics of the hadron-hadron interaction predict the universal 
character of the bp at the superhigh energies. 

3 Universality of cross sections at extremely large energies 

In the previous section we argued that rate of the increase of the size of the region where T(b, s) ~ 
1 should be followed by a rather steep drop of Y as given by single Pomeron exchange. Since the 
two scenarios of matching which we considered in the previous section give practically identical 
results for 10 3 < \fs < 10 6 GeV, we will consider here the soft matching dynamics which is 
obviously a more conservative way to estimate the approach to the unitarity limit. 

Imposing the condition of complete absorption at fixed b < bp and using the Pomeron ex- 
change formulae at larger b, we can build the partial amplitude so that it includes the continuity 
condition at the matching point 

T NN {h, s) = 9(bp - b) + exp [-(b - b F ) 2 /2B ]P \ 0(b - bp). (10) 
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Figure 1: Energy dependence of the black body regime cutoff parameter bp. Solid red line - 
estimate based on the soft matching approximation. Dashed blue line - prediction of the hard 
matching approximation. Here and in the following figures log s stands for logi (s/l GeV 2 ). 
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Since bp is universal at superhigh energies the only dependence on colliding particles is con- 
tained in the scale l/2Bjp for the impact parameter distribution. 

With the amplitude in Eq[lO| one can calculate the total cross section of the hadron-hadron 
interaction: 

a tot = 2 [ r{b,s)db = 2ir{b 2 F + 2B P ). (11) 



The slope of the t dependence of the elastic amplitude at t = is given by the formula: 

B = -J T(h - S)bdh= 2(V F + 2B F ) ' <12) 

So, 

-E- = 2(l ^ ) (13) 

AnB { (bl + 2Bj P y + 4Bl, ) - KV 

At accelerator energies where 2Bjp ^> bp we obtain relation: o /AnB « 1 . In contrast, at super 
high energies where bp ^> 2Bjp we obtain: 

a/AirB « 2(1 - AB 2 F /b F ). (14) 

Thus, at superhigh energies where 2a'ln(s/so) 3> £>o,iP the cross section and the slope of the t 
dependence at t = become the same for all colliding particles. The memory of the nature of 
colliding particles is lost as a consequence of the blackening of the interaction and the ratio of total 
cross sections for any colliding particles should be equal to one: 

crtot(hi, h 2 ) / crtotipp) -» !• (15) 

The same universality of structure functions is expected for superhigh energies because the 
interaction and essential impact parameters are increasing with energy and as the consequence of 
the U matrix unitarity condition [33 1: 

F 2A (x,Q 2 )/F m (x,Q 2 )^l. (16) 
A similar prediction holds for the total cross section of a photon-nucleus scattering: 

(Tja/vjp — » 1- (17) 



4 Glauber model cross sections for pA and AA collisions 

To visualize physical phenomena related to increase of the radius of a nucleon-nucleon interaction 
with energy and to learn whether the trend toward universality can be observed at the energies 
achievable at the accelerators or in cosmic ray we need to parameterize the energy dependence of 
the elementary pp cross section and the slope of t dependence for a wide range of energies. Up 
to the Tevatron energies (st ~ 4 • 10 6 GeV 2 ) we have chosen the parameterization of the total pp 
cross sections in the form satisfying the Froissart theorem : 

vtotipp) = *o(l + eln(s/s ) + e 2 /2\n(s/ s f) (18) 
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with e = 0.0808 . In the energy range where data are available this form is almost identical to the 
one suggested by Donnachie and Landshoff. The slope parameter B pp (t = 0) in this energy region 
is given by Pomeron exchange, 



Bpp(t = 0) = Bp = B 0t]P + 0.5 ln(s/l GeV 2 



(19) 



For higher energies we use as input the profile function T{b,s) which we built in the previous 
section in order to calculate the total, elastic and inelastic cross sections and the slope parameter of 
the iVW amplitude. The matching of the cross section and B NN is determined by using s = St as 
the reference point. The energy dependence of the elementary NN cross section and of the slope 
parameter Bnn(s) is shown in Fig|2| 

Now we can evaluate the total and absorption cross section for the proton-nucleus interaction 
at high energies in the conventional Gribov-Glauber model: 



oo 

23? J 1 - [l - J p(z,r t )T NN (b-r t ,s)dzdr t ]' 



db. 



(20) 



and, 



2a 



pp 



a 



pp 

tot 



p(z, r t )T NN (b - r t , s)dz dr t 



db. 



(21) 



Here p(r) = A~ 1 pa(y) is the single nucleon nuclear density normalized by the condition J p(r)dr = 
1. We calculated the nuclear density p& within the Hartree-Fock model with the effective Skyrme 
nucleon-nucleon interaction. Note that at intermediate energies this nuclear model provides a rea- 
sonable description of elastic proton and electron scattering off nuclei along the periodical table 
as well as the quasifree knockout of a nucleon in (e,e'p) and (p,2p) reactions without free param- 
eters J34|. In the above formulae we neglected correlations between nucleons, single and double 
inelastic diffraction. This is legitimate because for central collisions, where these approximations 
look suspicious, the cross section is close to black limit and therefore independent on details of the 
model. For the peripheral collisions these effects are small. 

The result of our calculations is shown in Fig|3]as a ratio of the total proton-nucleus cross sec- 
tion to the total cross section of the pp interaction as a function of the invariant energy s for p 16 
and p 208 Pb collisions. We present the ratio calculated in the approach of soft dynamics blacken- 
ing(solid red line) and in the hard regime (dashed blue line). The range between the two curves can 
be treated as a measure of uncertainty of our approximation. As we already discussed, both hard 
and soft regimes of blackening give close results up to the energies of the GZK limit. At higher 
energies, the hard mechanism leads to a faster approach to the universality regime. For comparison 
we also show the ratio found in the model neglecting the radius of the iViV interaction(dotted line). 
In the energy range where blackening is still a correction, neglecting the radius of the interaction 
leads to a stronger decrease of the (Jtrt/ a tot rat i° m me case of light nuclei. This indicates a more 
important role of peripheral interactions in the case of light nuclei. Even so, the universal asymp- 
totic is reached for light nuclei at extremely high energies, see Fig|4j The calculated absorption 
cross sections are shown in Fig|5] 

Accounting for the energy dependence of the radius of hadron-hadron interaction in the en- 
ergy domain where the NN cross section becomes large and the radius of the interaction becomes 
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Figure 2: Energy dependence of the total cross section and the slope parameter of iViV interaction. 
Solid line - the cross section and the slope parameter as dictated by the soft dynamics matching. 
Dashed line -hard mechanism of the blackening of NN interaction. Dotted line presents the slope 
parameter Bp due to the pomeron exchange 
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Figure 3: The dependence of a tot (pA) / <J to t(pp) on energy. Solid line - the ratio calculated with 
NN amplitude dominated by soft dynamics, dashed line - hard mechanism of the blackening of 
interaction. Dotted line is the ratio calculated neglecting by the radius of interaction. 
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Figure 4: The ratio of the proton-oxygen cross section to the proton-proton one as a function 
of energy. Calculation in the Gribov-Glauber model demonstrating the onset of the universality 
regime. 
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Figure 5: The absorption cross section of pA interaction as a function of energy. Dashed line - 
calculation of cross section neglecting the radius of interaction. 
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comparable to the radius of the nucleus reveals new effects beyond those usually associated with 
Gribov-Glauber shadowing. Scattering from the nucleus edge and from the meson "halo" of a 
nucleus should lead to a decrease in the dependence of the cross section on atomic number as com- 
pared to nuclear shadowing effects. It is evident that at asymptotic energies, where the interaction 
is already black in the range of the impact parameters considerably exceeding the radius of nu- 
cleus, the dependence on the atomic number in the nucleon-nucleus collisions should completely 
disappear. 

We also estimated how the energy dependence of the iViV interaction radius in the high energy 
domain will affect the total and inelastic nucleus-nucleus cross sections. These quantities are used 
in the study of the relativistic heavy ion collisions aimed to discover the new extreme state of the 
nuclear matter - the Quark-Gluon plasma. The calculation of the nucleus-nucleus cross sections 
in the Glauber-Gribov model is a rather complicated problem. Instead, for the rough estimates of 
the effect we used the generalization of the formulae of Bradt and Peters ll25l for the total cross 
section of scattering of two heavy ions: 

at? = 2<*{Rf{s) + Rf{s) - cf. (22) 

We use here the parameter c = 0.8Fm as found by Bradt and Peters, and we calculate the cross 
section of heavy ion collisions using the energy dependent nuclear radius R e ff(s) determined from 
the calculated proton-nucleus total cross section, 



f(s) = yJo£(s)/2n. 



R 

The energy dependence of the effective nuclear radius for lead is shown in Fig|6j The energy 
dependence of the PbPb total cross section is shown in Fig|7] where we also present the cross 
section given by the Bradt and Peters formula with the static nuclear radius(dotted line). The static 
nuclear radius was found from the calculated HF nuclear density (if one uses the empirical formula 
for the nuclear radius R A = r A 2 ^ 3 , the HF density of lead requires the value of r = 1.156). For 
comparison, we also calculated the total nucleus-nucleus cross section using the simplified Glauber 
approach expression [12611271 



1 - expl -T AB (h) 



(23) 



where 



T AB (h) = J dhi J dh 2 T NN (h 1 -h 2 ) J dzip A (bi, Zi) J dz 2 p B (b - b 2 , z 2 ). (24) 



We find significant corrections to the cross section of heavy ion collisions calculated using the 
Bradt-Peters model with the static nuclear radii already at RHIC energies. We also compare in 
FigE] these estimates of the total cross section to the results of calculation within the simplified 
Glauber optical model approach. In the range of energies 10 4 < s < 10 8 the Bradt-Peters formula 
underestimates the cross section due to the assumption of the sharp edges of nuclei, hence neglect- 
ing the contribution of the interaction of surface nucleons. However, at energies, close to the GZK 
limit the Bradt-Peters cross sections with an account of the energy dependence of the interaction 
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Figure 6: Increase with energy of the effective nuclear radius in pA collisions. Dotted line - 
neglecting the radius of interaction. 
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Figure 7: Cross section of nucleus-nucleus collisions. Dotted line - the Bradt-Peters cross section 
with static radii of nuclei; solid and short-dashed lines - the Bradt-Peters cross sections with en- 
ergy dependent nuclear radii as provided by hard and soft dynamics mechanisms correspondingly; 
dash-dotted line - optical Glauber approach without accounting for the energy dependence of the 
radius of iViV interaction; long dashed line - optical Glauber model with accounting for the energy 
dependence of the interaction radius 
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Figure 8: The average number of wounded nucleons in pA collision as a function of energy. 
Dotted line - calculation neglecting the energy dependence of the interaction radius. 
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Figure 9: The partial inelastic cross section as a function of the number of wounded nucleons. The 
yellow(light) bar is the cross section calculated neglecting the interaction radius. 
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radius should be reasonable for the total cross section as well as for the absorption cross section 
which is important for the interpretation of cosmic ray data. 

Now we want to demonstrate how accounting for the energy dependent radius of the iViV inter- 
action will change the distribution over inelastic collisions with increasing energy. We calculated 
the average number of wounded nucleons in nucleon-nucleus collisions 



v 



A ■ cr™ N ( 



(25) 



as a function of energy. Here, the inelastic cross section is calculated using the expression 



<7pAK s ) 



£ 

n=l 



(26) 



where the partial cross sections are given by formula 



cr n (s 



(A -l)\n\ 



dh 



o% N {s)TQ>,8) 



l-a% N (s)T(b,s) 



(27) 



with the generalized nuclear width function 



T(b, 



a 



tot I 



dr t T NN (r t -b,s) / dzp{v t 



(28) 



The energy dependence of the average number of wounded nucleons calculated with and without 
taking into account the increase of the radius of the iVW interaction with energy is shown in 
FigOO The effect is still small, on the level of 10%, at collider energies. At the same time, at 
asymptotic energies where the interaction becomes black in a wide range of impact parameters 
< v > \ s ^oo = A as a result of the universality of the collision amplitudes. We also find that, 
though < v > is weakly modified for collider energies, the partial cross sections are affected much 
more strongly, see Fig|9| In our calculations we neglected corrections due to inelastic diffraction, 
by the -momentum conservation law which are violated within the eikonal approximation. With 
the increase in radius of the nucleon-nucleon interaction the basic assumptions leading to Gribov- 
Glauber approximation, like the neglect of the simultaneous interactions with two nucleons which 
are located at different impact parameters, are violated since the radius of the NN interaction 
becomes comparable to the internucleon distance. 

The impact of these effects on the interpretation of the RHIC and future LHC data requires 
separate analysis which is beyond this paper. 



5 Discussion and Summary 

We demonstrated that the cross section of NN collisions at impact parameters growing as ln(s/so) 
reaches the black limit due to the hard dynamics. As a result, we argue that all hadronic and nucleus 
cross sections become equal at ultra high energies. We analyzed the role played by this effect for 
the total and absorption cross sections of the proton-nucleus collisions in a wide energy range. The 
effects are of the order 10% + 20% and, hence, should be taken into account in the future analyses 
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of the precision LHC data. They may also be relevant for interpretation of the cosmic ray data 
near the GZK limit. Dominance of hard dynamics at small impact parameters both in the elastic 
amplitude and in the structure of final states in the inelastic interactions at small impact parameters 
ID together with onset of the universality of the cross sections can be considered as new signals 
for a presence of the phase transition in NN interactions at central impact parameters. 
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